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Abstract 

The paper presents the QCD description of the small x behaviour of parton distribu- 
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in qualitative agreement with deep inelastic experimental data. 
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1 Introduction 



The recent measurements of the deep-inelastic (DIS) structure function (SF) F2 by 
the HI [|l| and ZEUS collaborations open a new kinematical range to study proton 
structure. The new HERA data show the strong increase of F2 with decreasing x, that 
contradicts to experimental data of the NMC || and E665 collaboration 0] at lower Q 2 
(Q 2 ~ lGeV 2 , which are quite flat at x ~ 1CT 2 . These lower Q 2 data are in the good 
agreement with the standard Pomeron or with the Donnachie-Landshoff picture || where 
the Pomeron intercept a p = 1.08, is very close to standard one a p — 1. However, the 
HERA data extracted at the larger Q 2 (Q 2 > lOGeV 2 ) are fitted very well (see || 0) 
by parametrizations of parton distributions (PD) contained the supercritical (or Lipatov) 
Pomeron. The interpritation of the fast changing of the intercept in the Q 2 region between 
Q 2 = lGeV 2 and Q 2 = lOGeV 2 (see [§]) is yet absent. There are arguments (see 0) in 
favour of one intercept or a superposition of two different Pomeron trajectories, one having 
an intercept of 1.08 and the one of 1.5 (see discussions of this problem in [§]). 

The aim of this letter is a possible "solution" of this problem in the framework of 
Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation [|10j. We will seek the "so- 



lution" g of DGLAP equation imposing the Regge-type behaviour (we use the parton 
distributions (PD) multiplied by x and neglect the nonsinglet quark distribution at small 

x): 

f a (x,Q 2 )~x- s f a (x,Q 2 ), (a = q,g) (a p = 1 + 5) (1) 

where f a (x,Q 2 ) is nonsingular at x — > and f a (x,Q 2 ) ~ (1 — x) u at x — > 10. Of 
course, we understand that the Regge behaviour (TJ) contradicts to the well-known double- 
logarithmical solution: ~ exp yJ(p(Q 2 )ln(l/x), where 4>{Q 2 ) is a known Q 2 -dependent 
function^ . However, we think that it is not possible to glue the both Regge-type be- 
haviour: Donnachie-Landshoff picture at low Q 2 and Lipatov Pomeron at large Q 2 , having 
the double- logarithmical (non Regge-type) solution at immediate Q 2 . 

In the case of the large 5 values (i.e. x~ 5 ^> 1) similar investigations were already 
done (see an< ^ ^he results are well known (see |l2j] for the first two orders of the 



perturbation theory, |TJJ for the first three orders and [JTJJ containing a resummation of 
all orders, respectively): 



f a (x,t) M a {l + 6,t) 



f a (x,t ) M a (l + 5,t ) 



(2) 



where t = ln(Q 2 /A 2 ), to = t{Q 2 = Qq) and M a (l + 5, t) is the analytical continuation of 
the PD moments M a (n,t) = Jq 1 dxx n ~ 2 f a (x, t) to the noninteger value "n = 1 + 5". Note 

2 We use the termin "solution" because we will work in the leading twist approximation in the range 
of Q 2 : Q 2 > lGeV 2 , where the higher twist terms may give the sizeable contribution. Moreover, our 
"solution" is the Regge asymptotic (|l|) with unknown parameters rather then the solution of DGLAP 
equation. The parameters are found from the agreement of the r.h.s. and l.h.s. of the equation. 

3 Consideration of the more complicate behaviour in the form x~ s (ln(l / x)) b i2g(\/ 4>ln(l/x)) is given 
inp 

^Vlore correctly, <j> is Q 2 -dependent for the solution of DGLAP equations with the boundary condition: 
fa(x,Qn) = Const at x — > 0. In the case of the boundary condition: f a (x, Qg) ~ exp ^Jln{l/x), <fi is lost 
(sec Q) its Q 2 -dependence 
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that recently the fit of HERA data was done in J7| with the formula for PD f q (x, t) very 
close to (H) and a very good agreement (the \ 2 P er degree of freedom is 0.85) is found 
for 5 = 0.40 ± 0.03. There are also fits |l6j of the another group using equations which 
are similar to (||) in the LO approximation. 

In this article we expand these results to the range where 5 ~ (and Q 2 is not large) 
following to the observed earlier (see [12, [T3] and Appendix) methodP] to replace the Mellin 



convolution by a simple product. 

Consider DGLAP equations and apply the method from |T3| to the Mellin convolution 
in its r.h.s. (in contrast with standard case, we use below a(Q 2 ) = a s {Q 2 ) / (An)): 

ji/«0M) = ~ \ Yl %ii a i x ) ® fifat) [a,b) = (q,g) 

i=a,b 

l - J2 7 ai (a,l + S)fi(x,t) + Oix 1 - 5 ) ( lab (a,n) = a^J {n) + a 2 7 % (n) + ...),{3) 



-a,b 

-)2 / \2\ 



where t = ln(Q /A ). The 7 a b(a,x) are the spliting functions corresponding to the 
anomalous dimensions (AD) 7o &(at, n) = Jq 1 dxx n ~ 2A f a b(a, x). Here the functions 7 a fc(a, 1 + 
5) are the AD r y a i } (a,n) expanded from the integer argument "n" to the noninteger one 
"1 + 5". The functions 7 a b(a, 1 + 5) (marked below as AD, too) can be obtained from the 
functions 7 a b(a, 1 + 5) replacing the term 1/5 by the one 1/5: 

1 ' \(l- V (u,5)x & ) (4) 



5 5 5 

This replacement M) appeares very naturally in the consideration of the Mellin convo 



lution at x — > (see |L3| and Appendix) and preserves the smooth and nonsingular 



transition to the case 5 = 0, where 

1 1 

~ = In g(u) (5) 

5 x 

The concrete form of the functions <p(v, 5) and g(v) depends strongly on the type of 
the behaviour of the PD f a (x, Q 2 ) at x — > and in the case of the Regge regime (HD they 
are (see Appendix ): 

<p(y*S) = ^(t + i-g)^ and = + !)-*(!). ( 6 ) 

where T(u + 1) and + 1) are the Euler T- and ^/-functions, respectively. As it can be 
seen, there is a correlation with the PD behaviour at large x. 

If 5 is not small (i.e. x~ s » 1), we can replace 1/5 by 1/5 in the r.h.s. of Eq.(^) and 
obtain its solution in the form fl2|). The case of small 5 values will be considered lower. 

The new points in our investigations are as follows. Note that the Q 2 -evolution of 
M a (l + 5,t) contains the two: "+" and "-" components, i.e. M a (l + 5,t) = £ i=± M*(l + 
5,t), and in principle each component evolves separately and may have independent (and 
not equal) intercept. 

5 The used formula (Eq.(2) from H) coincides with (||) in the leading order (LO) approximation, if wc 
save only f g (x,Q 2 ) in the r.h.s. of (|3|) (or put ^ qq = and "/ qg = formally). Eq.(||) and Eq.([|) from 
[Q have some differences in the next-to- leading order (NLO), which are not very important because they 
are corrections to the a-correction. 

6 The method is based on the earlier results |l7tl 
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2 Leading order 



Consider DGLAP equations for the "+" and "— " parts (hereafter s = ln(lnt/lnto)): 
^jtM = -^-Uoc,l + S ± )ft(x,t) + 0{x 1 - 5 \ (7) 
where 



7± 



1 



(igg + lgq) ± V (igg ~ 7 M ) + ^Igglgg 



are the AD of the "±" components (see, for example, ||i~8 



The AD 7_ (a, 1 + <L.) does not contain the singular term (see [|T^, [11] and below) and 
f~(x,t) have the solution in the form: 



/a~0Mo) 



-<W, where d + = 



2/3o 



The AD 7+ (a, 1 + 5+) contains the singular term and fai x i t) have the solution similar 
only for x~ s+ >> 1: 



/« + 0M) 



e -d+(l+8 + )s iix S+ >> 1 



(9) 



Both intercepts 1 + 5+ and 1 + 5- are unknown and should be found, in principle, 
from the analysis of the experimental data. However there is the another way. From 
the small Q 2 (and small x) data of the NMC and E665 collaboration [|J we can 
conclude that the SF F 2 and hence the PD f a (x, Q 2 ) have flat asymptotics for x — > and 
Q 2 ~ (1 -T- 2)GeV 2 . Thus we know that the values of 5 + and <5_ is approximately zero at 
Q 2 ~ lGeK 2 . 

Consider Eqs.(^) with 5± = and with the boundary condition f a (x,Qo) = A a at 
Ql = lGeV 2 . For the "— " component we already have the solution: Eq.(§) with 5- = 
and cL(l) = 16//(27/5 ), where / is the number of the active quarks and are the 
coefficients in the a-expansion of QCD /3-function. For its "+" component Eq.(^) can be 



rewritten in the form (hereafter the index 1 + 5 will be omitted in the case 5 



ln(-)-^-5 + (s) + ^-ln{Ai 
x as as 



1 

Wo 



7+(M-) - £<>)) + 7h 



where 7+ and 7 + are the coefficients of the singular and regular parts at 5 
1+ {l+5): 

1 4/ 

7+(l + 5) = 7+- + 7+, 7+ = -24, 7+ 



0): 

(10) 
of AD 



22 + — 
27 



The solution of Eq.(10) is 



f+(x,t) = A + a x^e- 2 ^ 



'.ID 



where 



(L 



7+ 
2/3 



4 

3' * 



2Po 



7+ 



4 . . 101 

-g(u) H 

3^ V ; 81 
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Hereafter the symbol ~ marks the case / = 3. 

As it can be seen from ([0]) the flat form 5+ = of the "+" -component of PD is very 
nonstable from the (perturbative) viewpoint, because d(5 + )/ds ^ 0, and for Q 2 > Ql 
we have already the nonzero power of x (i.e. Pomeron intercept a p > 1). This is in 
agreement with the experimental data. Let us note that the power of x is positive for 
Q 2 < Ql that is in principle also supported by the NMC [[J data, but the use of this 
analysis to Q 2 < lGeV 2 is open question. 

Thus, we have the DGLAP equation solution for the "+" component at Q 2 is close 
to Qq = lGeV 2 , where Pomeron starts in its movement to the supercritical (or Lipatov 
EH- 0) regime and also for the large Q 2 , where Pomeron have Q 2 -independent intercept. 
In principle, the general solution of (|7|) should contain the smooth transition between 
these pictures but this solution is absent f\ We introduce some "critical" value of Q 2 , Ql, 
where the solution flj^) is replaced by the solution ([□]). The exact value of Ql may be 
obtained from a fit of experimental data. Thus, we have in the LO of the perturbation 
theory: 

fa(x,t) = f~(x,t) + f+(x,t), f~(x,t) = A~ exp(-d_s) 

f+( ., j A+x^ex V (-d + s), HQ 2 <Q 2 C 

, ) \ f+(x,t c )e W (-d + (l + 6 c )(s-s c )), XQ 2 >Ql 1 ] 

where 

t c = t(Q 2 c ), s c = s(Q 2 c ), A- = A a - A + a and 

= (l-a)A q - aA g , A+ = aA g - eA q (13) 

and the values of the coefficients a, a and e may be found, for example, in |18 . 
Using the concrete AD values at 5 = and / = 3, we have 



1 44, + 9 A> A+ „ A , i A _ 4 ~ A 

27 Zn(I) - g(u) - ^ ' ^ " ^ + g Ar 243 ^ _ g{u) _ ^ 



Thus, the value of the "+" component of the quark PD is suppressed logarithmically 
and this is in qualitative agreement with the HERA parametrizations of SF F 2 (see 
211, |1) (in the LO F 2 (x, Q 2 ) = (2/9)f q (x, Q 2 ) for / = 3), where the magnitude connected 



with the factor x is 5 -r 10% from the flat (for x — > 0) magnitude. 



3 Next-to-leading order 

By analogy with the previous section and knowing the NLO Q 2 -dependence of PD 
moments, we obtain the following equations for the NLO Q 2 -evolution of the both: "+" 
and "— " PD components (hereafter s = ln(a(Ql)/a(Q 2 )),p = a(Ql) — a(Q 2 )): 

f a (x,t) = f~(x,t) + f+(x,t), f-(x,t)= A~ exp (-dS - d a __p) 

f+(r t) = / Atx ii+S+i + +p) exp (-d + ~s - d a ++ p), if Q 2 < Q 2 

Ja[ ' ) { f+(x,t c )ex V (-d + (l + 8 c )(s-s c )-d% + (l + 5 c )(p-p c )), if Q 2 > Q 2 C [ ' 



7 The form expf-S7+(l + 5)/(2/? )) coincides with both solutions: Eq.(g) if x d + >> 1 and Eq.([fl] 
when S = but it is not the solution of DGLAP equation. 



4 



where 



At 



s(Ql), Pc = p(Ql), "o = at(Ql), a c 



(l - a K£)At + a K a ± A^ 



d a = 
and K q ± 



d a ++ (H-) - g(y)) + d" 



7± T 



2/3o + 7± - 7t 



7±± 
7 , - 7 (o) 

± (0) 

7t - Jm 



7±A ~« 
K ± 



2$ 



(16) 



The NLO AD of the "±" components are connected with the NLO AD 7 ^. The 
corresponding formulae can be found in [TE| . 

Using the concrete values of the LO and NLO AD at 5 = and / = 3, we obtain 
the following values for the NLO components from fll5|) ,fll6|) (note that we keep only the 
terms ~ 0(1) in the NLO terms) 



d q 



16 



2C(3) + 9C(2) 



779 



108 



2800 - q 

-gT' d+ ~ 

"++ — gj ' "+- 



32 



28 

1.97, d 9 __ ~ + — ss 2.32 

81 

263 w „ s 372607i 



C(3) + -C(2) 
- + ^ - 12C(2) 



69984 
« -52.26 



-67.82 



(17) 



and 



A: 



20 



-Qf 



A? + g A 1 



3 

and A~ = A n - A 



1 44,(1 - 7.68a ) + 9A g (l - 8.72a ) 



9 



80 x 4 

Y a °) + 9 



27 

A, 



Hi) 



9A g - A q 



85 
108 



27/ 7l (I)_ e ( I/ ) 



85 
108 







692 



«o) 



It is useful to change in Eqs . (|T5|) - (p^8|) from the quark PD to the SF F 2 (x, Q 2 ), which 
is connected in NLO approximation with the PD in the following way (see fllfl): 

F 2 (x,Q 2 ) = (l + a(Q 2 )B q (l + 5))5 2 f q (x,Q 2 ) + a(Q 2 )B g (l + 6)5 2 f g (x, Q 2 ), (19) 

where S 2 = J2i=i If =< e / > is the average charge square of the active quarks: 5 2 = (2/9 
and 5/18) for / = (3 and 4), respectively The NLO corrections lead to the appearence 
in the r.h.s. of Eqs. ([15]) of the additional terms (l + aB±j + ctoB±j and the necessity 

to transform A^ to C^ 1 = F 2 ± (x, Q 2 ) into the input parts. The final results for F 2 (x, Q 2 ) 
are in the form: 



F 2 (x,t) = F 2 (x,t) + F+(x,t) 

F 2 ~(x,t) = C~ exp( y -d_5-d q __p)(l + aB-)/(l + a B-) 

c + x (d + s+d q ++ p) exp _ d q ++ p)(l + aB + )/(l + a B+ 

F+(x,t) = I F 2 + (x,t c )exp(-d + (l + 5 c )(s - s c ) - d q ++ {l + 6 c )(p - p c 
\ + aB+(l + 5 c ))/(l + a c B+{l + 5 C 



if Q 2 < Ql 

(20) 

if Q 2 > Ql 



where 



B 



± 



B n 



7± 

(o) "a-> 
Tg 9 



B g , C± = + 



with the substitution of A q by C = F 2 (x, Qq) into Eq.(|18D according 

C = (l+a B^5 2 A q + a B g 5 2 A 
For the gluon PD the situation is more simple: in Eq. 



(21) 

it is necessary to replace 
A q by C according 

For the concrete values of the LO and NLO AD at 5 = and / = 3, we have for 
(Revolution of F 2 (x, Q 2 ) and the gluon PD: 

F 2 (x,t) = F 2 ~{x,t) + F+(x,t), f g {x,t) = f~{x,t) + f+(x,t) 



F 2 -(x,t) 
F+(x,t) 

fg{x,t) ■ 

where 
C~ -- 

C + r 

At i 



32 8 
C~ exp(--s-1.97p)(l--a)/(l 

' C+x^ 



4 - | 2800 
3 s " 1 " 81 



p>«Kp(-f (<?(,/) + + 



2800 
81 



el* 7 . 



+ 67.82)p) 



6[Jn(±) - - |]a)/(l + 6[in(J) " - ifQ 2 < 



F 2 + (x,t c ) exp(-d+(l + 5 C )(5 
1 + a£+(l + 4))/(l + a c 5+(l + 8 C ) 
32 



^c)(p-p c 



22) 



ifQ 2 > Q\ 



A- exp(-— s-2.32p) 

A^-P^-Wp) exp (-f (g(z/) + 1)5 + C-^g(u) + 52.26)p), 
f+(x,t c )exp + 5 c )(s - s c 



4+(l + <y c )(p-p t 



ifg 2 <Q2 

_ ~ , 23 

ifg 2 >gf ; 



c 

2 

27 



26ao 



A + 2C 



^9 _ A g 



+ 



and 

10.50a ) + 2C(1 - 8.55a ) 



A, 1 



9 



■a 



2C 1 



Hi) - e{v) - n ; 

(><)2 , 2 2A g (l - f^ ) - C(l - f^ ) 



51 



■a 



27 



(24) 



Let us give some conclusions following from Eqs.(|22|)-(]24]). It is clearly seen that 
the NLO corrections reduce the LO contributions. Indeed, the value of the supercritical 
Pomeron intercept, which increases as ln(ao/a) in the LO, obtains the additional term 
~ (qjo — a) with the large (and opposite in sign to the LO term) numerical coefficient. 
Note that this coefficient is different for the quark and gluon PD and this is in agreement 
with the recent MRS{G) fit in and the data analysis by ZEUS group (see f22]| ). The 
intercept of the gluon PD is larger then the quark PD one (see also || j2~2fl). However, the 
effective reduction of the quark PD is smaller (which is in agreement with W.-K. Tung 
analysis in [23]), because the quark PD part increasing at small x obtains the additional 
(~ cto but not ~ 1/lnx) term, which is important at very small x. 

Note that there is the fourth quark threshold at Q 2 h ~ lOGeV 2 and the Q\ h value may 
be larger or smaller to Q 2 one. Then, either the solution in the r.h.s. of Eqs. (|20|j22]j23|) 
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before the critical point Q 2 and the one for Q 2 > Q 2 C contain the threshold transition, 
where the values of all variables are changed from ones at / = 3 to ones at / = 4. The 
a(Q 2 ) is smooth because — > (see also the recent experimental test of the flavour 
independence of strong interactions into p4||). 

For simplicity here we suppose that Qf h = Q 2 and all changes initiated by threshold 
are done authomatically: the first (at Q 2 < Q 2 ) solutions contain / = 3 and second (at 
Q 2 > Q 2 ) ones have / = 4, respectively. For the "— " component we should use Q 2 h = Q 2 C1 
too. 

Note only that the Pomeron intercept a p = 1 — (d + s + d q ++ p) increases at Q 2 = Qf h , 
because that agrees 

_ 1 / IKQLQD - ^rP(Qth,Qo), HQ 2 <Q 2 C 

P ~ ~ 1 l.U5(Q 2 th ,Q 2 ) - 38.11p(Q 2 th ,Q 2 ), if Q 2 > Q\ 

with results obtained in the framework of dual parton model. The difference 

Aa p = 0.irs(Q 2 th ,Ql)-3.55p(Q 2 h ,Q 2 Q ) 

dependes from the values of Q 2 th and Q\. For Q 2 h = lOGeV 2 and Ql = lGeV 2 it is very 
small: 

Aa p = 0.012 



4 Discussion 

Let us summarize the obtained results. We have got the DGLAP equation "solution" 
having the Regge form ([]]) for the two cases: at small Q 2 (Q 2 ~ lGeV 2 ), where SF and 
PD have the flat behaviour at small x, and at large Q 2 , where SF F 2 (x, Q 2 ) fastly increases 
when x — > 0. The behaviour in the flat case is nonstable with the perturbative viewpoint 
because it leads to the production of the supercritical value of Pomeron intercept at larger 
Q 2 and the its increase (like 4/3 ln(a(Ql) / a(Q 2 ) in LO) when the Q 2 value increases^. The 
solution in the Lipatov Pomeron case corresponds to the well-known results (see [0, [14|. (7j) 



with Q 2 -independent Pomeron intercept. The general "solution" should contains the 
smooth transition between these pictures. Unfortunately, it is impossible to obtain it 
in the case of the simple approximation ([[]), because the r.h.s. of DGLAP equation (0) 
contains both: ~ x~ s and ~ Const terms. As a result, we used the two above "solutions" 
gluing them at some point Q 2 . 

Note that our "solution" is some generation (or an application) of the solution of 
the DGLAP equations in the momentum space. The last one has two: "+" and "— " 
components. Our conclusions are related to the "+" component, which exhibits the basic 
Regge asymptotic behaviour. The Pomeron intercept corresponding to "— " component, is 
Q 2 -independent and this component is the subasymptotical one at large Q 2 . However, the 
magnitude of the "+" is suppressed like l/ln(l/x) and a(Ql), and the subasymptotical 
"— " component may be important. Indeed, it is observed experimentally (see El], 
Note, however, that the suppression ~ «(<5o) ^ s really very slight if we choose a small 
value of Qq. 

8 The Pomeron intercept value increasing with Q 2 was obtained also in pH, [27[|. 
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Our "solution" in the form of Eqs.fl2"2|)-(P1D is in very good agreement with the recent 
MRS(G) fit || and with the results of at Q 2 = 15GeV 2 . As it can be seen from 
Eqs.(|22|),(f23|), in our formulae there is the dependence on the PD behaviour at large x. 
Following PS| we choose v = 5 that agrees in the gluon case with the quark counting rule 
29]. This v value is also close to the values obtained by CCFR group [30] iy = 4) and 
in the last MRS{G) analysis || (v = 6). Note that this dependence is strongly reduced 
for the gluon PD in the form 



f g (x,Ql) = A g (v)(l-x) 1 



if we suppose that the proton's momentum is carried by the gluon, is ^-independent. We 
used A g {5) = 2.1 and F 2 (x, Q 2 ) = 0.3 when x -> 0. 

For the quark PD the choice v = 3 is more preferable, however the use of two different 
v values complicates the analysis. Because the quark contribution to the "+" component 
is not large, we put v = 5 to both: quark and gluon cases. Note also that the variable 
v(Q 2 ) has (see j^TJ) the Q 2 - dependence determined by the LO AD 7^. However this 
Q 2 -dependence is proportional to s and it is not important in our analysis. 

Starting from Ql = lGeV 2 (by analogy with p2| ) and from Ql = 2GeV 2 , and using 
two values of the QCD parameter A: the more standard one (Ag = 200 MeV) and 
(Ag= 255 MeV) obtained in [H], we have the following values of the quark and gluon 



PD "intercepts" S a 
if Q 2 = 1 GeV 2 



(d + s + d a ++ p) (here A|^| is marked as A): 



if Ql 



Q 2 


W) 

A = 200MeT/ 


A = 200MeV 


A = 255M e y 


W 2 ) 

A = 255MeV 


4 


0.191 


0.389 


0.165 


0.447 


10 


0.318 


0.583 


0.295 


0.659 


15 


0.367 


0.652 


0.345 


0.734 


2 GeV 2 


Q 2 


W) 

A = 200MeV 


W) 

A = 200M e y 


A = 255M e y 


A = 255^^ 


4 


0.099 


0.175 


0.097 


0.198 


10 


0.226 


0.368 


0.227 


0.410 


15 


0.275 


0.438 


0.278 


0.486 



These values of 5 a are above at 4 GeV 2 those from M. Because we have the second 
(subasymptotical) part, the our effective "intercepts" have smaller values. 

Note that as input we can use the Eq.(|l]) with 5 = e = 0.08, which corresponds to 
Donnachie-Landshoff value |33[ of the Pomeron intercept. For this purpose we should 
represent the value \je (see Eq.(f|)) as the series (here \1/ [y) = 4-^(u)): 



In- 



r(z/ + i)r(i-e) 

V{y + \-e) 
^(u + l) - 



In— 

x 



1) - *(1) 



+ 



(25) 
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and save only the first term in the r.h.s. This is possible if the second term in the r.h.s. 
of Eq.fl2~5|) is negligibly small, i.e. 



In— 

x 



, i 

< ln- 

x 



1) - *(1) 



or for v — 5 



1.5 



10" 12 < x <€. 9.6 ■ 10" 2 



This x range corresponds to the range considered here. The right boundary slowly 
depends on the S value: it is determined by the large x SF behaviour, i.e. by the v value. 
The left boundary strongly increases with increasing of the 5 value, that leads to the 
impossibility to apply the expansion (HE) for the large 5 values. 

In the case 5 = e all our conclusions are not changed except the relation a p = 1+5(Q 2 ), 
which transforms to the form a p = 1 + e + S(Q 2 ), i.e. the value of the Pomeron intercept 
slightly increases. 

In support of our analysis we represent the Fig. 3 of paper || and add it by our 
results for 5 q (Q 2 ) of above table at Q 2 , — 1 GeV 2 . The dashed-dotted curve represents 
the intercept of the Pomeron trajectory a p (Q 2 ) which was obtained in |2(| as the result 
of a fit of experimental data. This analysis was done in the framework of Regge-type 
behaviour of DIS SF, that corresponds to the start point of our consideration, too. The 
authors of |2Bj obtained very strong Q 2 -dependence of the Pomeron intercept in the region 



of 1 GeV 2 < Q 2 < 10 GeV 2 and approximate Q 2 -independent its values a p ~ 1.05 and 
a p ~ 1.4 at Q 2 < 1 GeV 2 and Q 2 > 10 GeV 2 , respectively. The solid curves represent 
our values of Q 2 -dependence of Pomeron intercept 



a P (Q 2 



1.05 + 5 q (Q 2 ) 
1.05 + 5 q (Q 2 c ) 



if Q 2 < Ql 
if Q 2 >Ql 



(26) 



where we choose a p (Ql) = 1.05 and Ql = 15GeV 2 . As can be seen in the Figure both 
the results are in very good agreement. 

As a conclusion, we note that BFKL equation (and thus the value of Lipatov Pomeron 
intercept) was obtained in |l!| in the framework of perturbative QCD. The large-Q 2 
HERA experimental data are in the good agreement with Lipatov's trajectory and thus 
with perturbative QCD. The small Q 2 data agrees with the standard Pomeron intercept 
a p = 1 or with Donnachie-Landshoff picture: a p = 1.08. Perhaps, this range requires 
already the knowledge of nonperturbative QCD dynamics and perturbative solutions (in- 
cluding BFKL one) should be not applied here directly and should be corrected by some 
nonperturbative contributions (see 



5 Conclusions 

Thus, in our analysis Eq.([l|) can be considered as the nonperturbative (Regge-type) 
input at Ql ~ lGeV 2 . Above Ql the PD behaviour obeys DGLAP equations. Under 
the action of perturbative QCD the Pomeron splits into two components. The intercept 
of the "— " component is Q 2 independent and this component is the subasymptotical 
one at larger Q 2 values. The Pomeron corresponding to the "+" component moves to 
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the supercritical regime and tends to its perturbative value. Above some Q 2 C1 where its 
perturbative value is already attained, the Pomeron intercept keeps a constant value. Our 
analysis supports the idea |2(| |27|] about the one effective Pomeron having a Q 2 -dependent 
intercept, however the character of the Q 2 -dependence is different. 

The application of this approach to analyse small x data and taking into account a 
resummation of all c^-orders of perturbative QCD (by analogy with |T5[) invite further 
investigation and will be considered in future. 
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6 Appendix 



Here we present the illustration of the method to replace the convolution by simple 
product at small x. More detailed analysis can be found in 
1. Consider the basic integral 

J s (a, x)=x a * ip{x) = f 1 ^ y a <p{-), 

J* y y 

where <f(x) = Ax~ 5 (l — x) u = x~ 5 (p(x). Expanding (p{x) near <p(0) , we have 



J s (a,x) = x- d I dyy^ 5 - 1 
1 



^(o) + -fp«(o) + ... + i(-) ^ (fc) (o) + 



x 



x 



a + 5 
1 



0(O)+O(x) 



(Al) 



a + 5 



a + S - 1 



^ (1) (o) + ...+ 



k\ a + S — k 



v (k \o) + ... 



The second term on the r.h.s. of eq.flAI]) can be summed, and J$(a } x) has the following 
form 



Js(a, x) = x 5 
tw< 

a) a > 1 



a + 6 

Consider two important cases: 



<f(0) + O{x) 



1 (1 + v — a — o) 



Js(a, x) = x~ 



<p(x) + 0(x 



2-<5^ 



b) a = 

J 5 (0,x) = x~ s 



1 



<p(0) + O(x) 



+ 



a + 5 

r(-<j)r(i + i/ 



i 



^(0)=x- 5 T ^(x) + O(x 1 - 5 ) 
1 (1 + v — o) S 
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where 



1 



i> + i)r(i-£) s 
r(v + i-8) 



2. Consider the integral 

I s (x) = K(x) * <p(x) = 

and define the moments of the kernel K{y) in the following form 



* y y 



dy y n - 2 K(y) 



In analogy with subsection 1 we have 

h{x) = x~ s f 1 dy y 5 " 1 K{y) 0(0) + - 0«(O) + . . . + ± (-) ' 
Jx y k\ \y j 

= x- s [K 1+5 <p(0) + O(x)\ 



(A2) 



x 



N 1+5 (x) <p(0) + N 5 (x) ^W(O) + . . . + - N 1+S _ k (x) £<*>(0) + ... 



where 



N v (x) 



dy y v 2 K(xy) 



The case K 1+s = l/(a + S) corresponds to K(y) = y a and has been already considered 
in subsection 1. In the more general cases (for example, K 1+ g = ^(8) +7) we can represent 
the "moment" K 1+ s as the combination of singular and regular (for <5 ^ 0) parts, i.e. 
Ki + s = —1/8 + \1/ (1 + 8) + 7. For the singular term the analysis from subsection 1 may be 
repeated. As the regular part can be represented by series of the sort J2 m =i l/(a + 8 + m), 
then any additional contributions from term Ni + s(x)(p(0) to any term of the series, are 
not necessary. 

So, for the initial integral at small x we get the simple equation: 

h(x)= x- 5 K l+S ^(x) + 0(x 1 - 5 ) 
where Ki+s coincides with K± + s after the replacement 1/8 — > 1/8. 
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Figure 1. The intercept of the Pomeron trajectory a p (Q 2 ) (dashed-dotted line) as 
obtained from the ALLM parametrization (see || [2(J). The dotted lines show the uncer- 
tanty of the fit. The solid curves represent the values of a p {Q 2 ) from Eq.([26l). 
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